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Abstract We present a class of conformally flat solutions of the Einstein’s field 
equations for spherical systems undergoing gravitational collapse accompanied 
with radial heat flux. The interior space-time of the collapsing matter is chosen 
to be conformal to the Robertson-Walker space-time so that its Weyl tensor van¬ 
ishes identically. The conditions which ensure the smooth matching of the interior 
space-time of the collapsing system, across its boundary, with the exterior space- 
time of a radiating star which is described by Vaidya metric lead to an equation 
governing its overall subsequent evolution. New solutions of this equation have 
been shown to provide physically viable models of collapsing stars. 

Keywords General relativity • Conformally flat space-time • Gravitational 
collapse • Heat flux. 


1 Introduction 

Study of relativistic models of gravitationally collapsing systems has been receiv¬ 
ing considerable attention since Oppenheimer and Snyder pQ first examined the 
collapse of a homogeneous spherical distribution of matter in the form of dust. 
It is known that a star which loses its equilibrium, collapses under self gravity, 
ultimately culminating in the state of space-time singularity. The state of singu¬ 
larity may be a naked one, in the sense that light rays emanating from it can reach 
any far-away observer or a black hole implying that the light rays from it may be 
trapped within a surface called its event horizon. Our understanding about the 
final outcome of any collapsing stellar configuration is known to rely on a couple 
of adopted conjectures only in this regard |21I51 - Accordingly, studies of realistic 
models of gravitationally collapsing systems are expected to play a key role in 
promoting our understanding of the nature and dynamics of the collapsing stellar 
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systems. Many studies aimed at examining various aspects of gravitationally col¬ 
lapsing stellar systems of different kinds of matter distributions have been carried 
out on varying types of space-time background. The space-time in the exterior re¬ 
gion of a radiating fluid sphere is known to be described by Vaidya[3] metric. The 
formulation of the boundary conditions smoothly joining the interior space-time- 
metric of the collapsing matter across its boundary with the appropriate form of 
Vaidya metric of the exterior space-time, as first proposed by Santos [5] has given 
a tremendous impetuous for studies in this direction. 

In the present work, we have proposed a relativistic model of a spherically 
symmetric matter source, whose collapse is accompanied with dissipation in the 
form of radial heat flux following Santos [5] ■ The background space-time metric is 
chosen so as to have vanishing Weyl tensor implying its conformal flatness. Confor¬ 
mally flat space-times, in the context of radiating fluid spheres, were first studied 
by Som and Santos[7j. Later, the most general class of conformally flat solutions 
for a shear-free radiating star were obtained and examined by Maiti[8], Modak[9], 
Banerjee et af[10], Patel and Tikckar|l 1 j. Schafer and Goennerjl2l and Ivanov[T5|. 
Herrera et al[ ll, have critically examined models of shear-free collapsing fluids 
accompanied with dissipation of heat on the space-time background subject to the 
constraint that the associated Weyl tensor should vanish. In the present approach, 
the geometry of the background space-time of the non-adiabatically collapsing 
matter is chosen to be conformal to that of the space-time of Robertson-Walker 
metric so that the Weyl tensor is known to vanish since the Robertson-Walker 
space-time is already conformally flat. 

The paper has been organized as follows. In Sec.0 the Einstein field equations 
for non-adiabatically collapsing spherical matter distributions with radial heat 
flow have been formulated on the conformally flat space-time background with 
pre-assigned metric form. The pressure isotropy relation in this approach is found 
to lead to two classes of solutions for the system. The junction conditions which 
smoothly join the interior space-time with the exterior Vaidya[I] metric across 
the boundary surface £ have been obtained in Sec. [3] In Sec. [H the explicit 
expressions for the physical quantities and the constraints on the model-parameters 
in view of physical requirements have been obtained. In Sec. 0 specific models 
have been examined by considering a couple of solutions of the equation governing 
the evolution of the collapse and their physical viability is shown using graphical 
methods. Sec. [G] contains a discussion accompanied with some concluding remarks. 


2 Conformally flat radiating collapse model 


We begin with the space-time describing the interior of a spherically symmetric 
collapsing matter fluid with radial heat flux with its metric 


ds = 


A 2 (r,t) 


dt 


dr 2 

1 — kr 2 


— r 2 {d0 2 + sin 2 Od<f> 2 ) 


( 1 ) 


The evolution with time of the configuration is governed by the function A(r, t) 
determined by relativistic field equations. The Weyl tensor for the space-time of 
metric m will vanish indicating its conformal flatness since the space-time is also 
conformal to the Robertson-Walker metric which follows on setting A(r,t) = a(t). 
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The matter content of the collapsing object in the presence of heat flux is 
described by the energy momentum tensor 

T.f = (p + p)uiu j - pS^ + qiu j + q j Ui, (2) 

where p is the energy density, p is the isotropic fluid pressure, u l = (A, 0,0,0) 
is the 4-velocity of the fluid and q l = (0, <7,0,0) is the radially directed heat flux 
vector. 

Subsecfuently, Einstein field equations relating metric potentials of m with the 
dynamical variables of © is a system consisting of the following four equations: 

p = 3 kA 2 + 3 A 2 - 3(1 - kr 2 )A' 2 + 2(1 - kr 2 )AA" + -(2 - 3 kr 2 )AA', (3) 

r 

p = -kA 2 + 2AA - 3 A 2 + 3(1 - kr 2 )A' 2 - -(1 - kr 2 )AA\ (4) 

r 

p = -kA 2 + 2AA - 3 A 2 + 3(1 - kr 2 )A' 2 - 2(1 - kr 2 )AA" 

--(1 - 2kr 2 )AA\ (5) 

q = —2(1 — kr 2 )A 2 A', (6) 


where an overhead prime ( ; ) and an overhead dot (.) denote differentiations with 
respect to r and t, respectively. 

Eqs. © and (0) lead to 


AT_ 

A 


A' 1 
rA (1 — kr 2 ) 


= 0 . 


(7) 


Eq. © usually referred to as pressure isotropy condition is easily integrable and 
admits two solutions: 

Case I: k = 0: 

A(r,t) = £,{t)r 2 (8) 

Case II: k ^ 0: 

"4( r , t) = £ (t) \/5T- kr 2 +r](t). (9) 

On setting £ = 0 one obtains the FRW model of the universe while on setting £ = 
a constant, one obtains a perfect fluid space-time solution without any dissipation 
or heat flow. Accordingly, solutions obtained on setting time dependence of the 
functions £(t) and p(t) appropriately will lead to models of collapsing fluid systems 
in the presence of heat flux. 

Note that, in the context of shear-free radiating collapse, the conformally flat 
space-time metric obtained earlier by Herrera et aim has the form (in coordinates 
(t,r,0,<f>)) 


ds 2 _ = 


[X(t)r 2 + l ] 2 d p 

[£(f> 2 + p(t)} 2 


1 

[£( t ) f 2 + f )(£)] 2 


[dr 2 + r 2 (d0 2 + sin 2 0dcj> 21 \ . 


( 10 ) 


Obviously, the solution © is a sub-class of the Herrera et aim model which can 
be obtained by setting x(f) = 0 in Eq. (1101) . Further, if we set X = k and make the 
following coordinate transformations 


r 

(1 + y/1 — kr 2 ) ’ 
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the solution d9j) can be obtained from m by identifying the parameters 77 = fj+^/k 
and £ = fj — £/fc. Hence, the solution d9j) also turns out to be a sub-case X = a 
constant of Eq. m- Since all conformally flat space-times are conformally related, 
this is expected. 

It should be stressed here that even though the classes of solutions obtained in 
this paper are special cases of the conformally flat solution of na, the solutions 
have been obtained by adopting a different approach. The conformally flat solution 
of [T2] has been obtained by equating the Weyl tensor of the associated space-time 
to zero while the solutions reported in this paper have been obtained by impos¬ 
ing the pressure isotropy condition in a space-time with geometry conformal to 
a conformally flat Robertson-Walker space-time. The space-time geometry of the 
non-adiabatically collapsing fluid sphere continues to be conformal to the homo¬ 
geneous geometry of the RW space-time as the collapse proceeds. In addition, the 
two classes of solutions have the following distinctive features. If £(t) = ( k/4)fj(t ), 
x(t) = (fc/4), on choosing a new time coordinate dT = ), one finds that the 

space-time of dini) corresponds to that of a FRW fluid sphere without heat flux. 
In our model, the space-time degenerates to that of a FRW fluid sphere without 
heat flux when £(f) = 0. It is noteworthy that our choice of curvature coordinates 
brings out to attention the role of the curvature parameter k of RW space-time 
explicitly in the evolution of the collapse. 


3 Junction conditions 

The space-time exterior to the collapsing matter source will be filled with radiation 
and is appropriately described by the Vaidyajlj metric for a radiating star 



( 11 ) 


The interior space-time 0 should smoothly be joined with that of the Vaidya 
metric m across the boundary surface S separating the interior (r < m) and the 
exterior (r > rs) space times. Following Santos[5], we shall stipulate the boundary 
conditions which ensure continuity of (i) the metrics of the interior and exterior 
space times with that on 27, and (ii) the extrinsic and intrinsic curvatures of S. 

Let gij be the intrinsic metric to U such that 


ds% = gijdC , z = 1,2,3 


( 12 ) 


and be the metric corresponding to the exterior (+)/interior (—) regions so 
that 


ds± = g^pdx±dx±, o = 0,1,2,3. (13) 


(13) 


Then the first condition implies 


(dsi.) s = (ds^s = ds%. 


(14) 


The second condition implies 


(K+) s = (K~ j )s, 


( 15 ) 
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where, 


K^ 


d 2 X% 


zh f-\d. 
fiis 


dx £ d x "± 

dE, 1 d& ' 


(16) 


In (1161) . n„ denote the components of the normal vector to S in x± coordinates. 
We express the intrinsic metric on S as 


ds% = dr 2 — R 2 (r)(d9 2 + sin 2 Odcj) 2 ). 


(17) 


Using the condition m for the intrinsic metric ED of the interior space-time i, 
we obtain 


f; =A(r s ,t ), (18) 

rs = R(T)A(rs,t). (19) 

In the interior space-time 

f(r,t) = r ~r s = 0. (20) 

The vector with components is orthogonal to U. Accordingly the unit normal 
to S in coordinates x~ is obtained as 

““ =1 °’aTT^SV 0 ’ 0 *' (21) 

The intrinsic curvatures corresponding to the interior region have the explicit 
expressions 


Kgg = 


A\ 1 - kr 2 ) 1 ' 2 

r r 2 A' 

A ~ IE 


J u 

(1 -kr 2 ) 1 / 2 


( 22 ) 

(23) 


Using the condition ED for the metric ED of the exterior space-time ED , we 
obtain 


r s (v) = R(t), 

( C =( 2 £ +i -^1 

For the exterior region, we have 


/( r,v) = r - r(u) = 0, 


and 


df 

dx+ 




(24) 

(25) 


(26) 

(27) 


The unit normal to S is 
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The expressions for the extrinsic curvatures corresponding to the exterior region 
are found to be 


K 


+ __ 

TT 


Kte = 


( d 2 v\ f dv\ 1 m dv 
dr 2 J \dr J r 2 dr 



(29) 

(30) 


where Eq. ([25]) and its derivative were used to derive Eq. ([29]) . Using Eqs. m , 
uni), (Hi and (1251) in Eqs. (1231) and o and imposing the condition 

= (RJ 9 ) s , (31) 


we obtain 


m(v) = 


2 A 


, A' 


kr + 2r(l — kr )~— + r [ — ) — r (1 — kr ) ( 


A 


A 


A' 


A 


■ (32) 


It follows from Eq. (1321) . that the mass within the sphere of radius r s at any 
instant t can be expressed as 


^ s r 
m(r ’ <) = 2 A 


kr 2 + 2r(l — fcr 2 )^- + r 2 


— r 2 (l — kr 2 ) 


■ (33) 


Using Eqs. ]5l), ]6j), (fl8l) . (fl9l) . (l22l) . (l23l) . (l24l) . (IM1) . (l30l) and (l32l) and imposing 
the condition 

(K+) s = (K- r ) s , (34) 


we finally obtain 


s 
P = 


A\/l — kr 2 _ s 


(35) 


4 Behaviour of the physical quantities 

We shall consider the cases k = 0 and k 0 separately. 


4.1 Case I: k = 0 

The models of collapsing fluid configurations of this class are characterized by the 
following expressions for energy-density, pressure and heat flux, respectively: 

P = 3 [ 4^77 + (r 2 i + t)) 2 ] , (36) 

p = 12r 2 £ 2 - 8£(r 2 £ + rj) - 3(r 2 £ + rff + 2(r 2 £ + p)(r 2 £ + rj), (37) 
q = -4r(r 2 £ + rjfi. (38) 

The inhomogeneities in matter density distribution are observed to be linked with 
the dissipative processes and disappear when the heat flux vanishes. 
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The total mass of fluid contained within the spherical region of radius rs at 
any instant t has the explicit expression 


m(r, t) 


s r 3 [4g?7 + ( r 2 £ + ?)) 2 ] 
2 (r 2 £ + 7/) 3 


(39) 


The parameter of expansion is 

0 = uf a j = — 2>A = —3 (r 2 £ + 77 ). 

In view of Eqs. (1371) and (1381) the condition (13511 . is equivalent to 

2 {r%£ + rf)(r%£ + rf) - 3 (r%£ + r/f + 4 rz{r%£ T ??)£ 
-8 + 77) + 12rk 2 = 0. 


(40) 


(41) 


This relation governs evolution of the physical quantities. Eq. m is a highly non¬ 
linear equation relating two unknown arbitrary functions (£(f), 77 (f)). A suitable 
choice of one of the two functions is essential to determine the other. The functions 
are further constrained by the requirements of (weak) energy conditions ( p , p, q > 
0) and other physical requirements. For a collapsing sphere 0 should be negative 
implying (r 2 f + 77 ) >0. Positive heat flux is ensured if £ < 0. From (1361) . we 
conclude that dp/dr = 12r£(r 2 £ + 77 ). Since, (r 2 ^ + if) > 0 and £ < 0, it follows 
that dp/dr < 0. The implications of other physical requirements can be examined 
only by adopting numerical procedures. 


4.2 Case II: k ^ 0 


The energy-density, pressure and heat flux of the collapsing fluid configurations of 
this class have the following expressions: 


P = 3 


—kf 2 + krj 2 + (1 — kr 2 )f 2 + 2\J\ — kr 2 £i] + 77 2 , 


p = [3fc 2 r 2 £ 2 — k(y /1 — kr 2 £ + 77 ) 2 + 4fc£((l — kr 2 )f + \J 1 — kr 2 rf) — 
3(\/l — kr 2 ^ + 1)) 2 + 2(\/l — kr 2 £ + 77 ) (\/1 — kr 2 £ + 77 )], 
q = 2 kr\J 1 — fcr 2 (\/ 1 — &r 2 £ T 77) 2 £. 


(42) 

(43) 

(44) 


The presence of dissipation in this case also is observed to give rise to inhomo¬ 
geneities in matter density distribution of the collapsing fluid. 

The total mass of the fluid enclosed within the spherical region of radius rs of 
the configuration of this class at any instant t is 


2(Vl - A:r 2 e + 77)3 
The expansion parameter is 


fc(-f 2 + rf 2 ) + (1 - kr 2 )fp + 2\J 1 - kr 2 f,i] + 77 2 . 

(45) 


0 = —3(\/l — kr 2 £ + if). 


(46) 
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In view of Eqs. m and (14411 , the condition (13511 . is equivalent to 

2 (\A “ kr st + 7 ?)(\/ 1 “ + i i)~ 2 fer^(^/l - fcr|,£ + 77)^ 

+3 r%k 2 f - 3(^/l - kr%i + t)) 2 - k(^j 1 - kr\i + r;) 2 

+4£fc[(l - + 0. - fcr|?7] = 0, (47) 

which is the relation governing evolution of the various physical quantities. Again, 
one needs to make a suitable choice of one of the two arbitrary functions to deter¬ 
mine the other. The functions are further required to comply with the requirements 
of energy conditions ( p , p, q > 0) and physical plausibility requirements. In this 
case, O will be negative if \/l — kr 2 £ + 77 > 0. Positivity of heat flux is ensured 
if £ > 0 for k > 0 and f < 0 for k < 0. The implications of other physical 
requirements can be examined using numerical procedures. 


5 Particular models 

In view of the complexity of the relation governing the evolution with time of the 
collapse, we have used two approaches for examining the physical plausibility of 
the evolution. 


5.1 Approach: I 


The evolution of the collapsing systems is governed by Eqs. m and (|47|) . In order 
to analyze physical viability of the model, it is essential to solve these equations. 
Note that the metric HD is conformal to the Roberson-Walker metric for which 
we have k = 0, ±1. For k = +1, we have 0 < rs < 1. However, there is no such 
restriction on rs ( 0 ,oo) for k = 0 and k = — 1 . 

We consider a particular case k = +1 for which the maximum value of the 
coordinate parameter is r^j = 1. Substituting these values in Eq. (HU), we obtain 

27777 - 2rji - 3 77 2 + 3 £ 2 - ? 7 2 = 0. (48) 


We need to specify one of the two unknowns to determine the other. Accordingly, 
we assume 



(49) 


where, n is a constant. Solution of Eq. <135)1 is then obtained as 


£ 


e V3t(^ + 


3t 2 (\/3e^3‘ + 3 Q) 


(50) 


where, Q is an integration constant. 

Thus, a particular model has been developed which can be utilized to analyze 
evolution of the physical parameters. For this particular model, Fig. |T] and [2] indi¬ 
cate evolutions of the surface density and heat flux, respectively, beginning with 
zero values at t = —00 when r^j = 1. Fig. [3] indicates that the mass decreases 
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monotonically as t —> 0. In Fig. [4] the expansion rate (<9) has been plotted which 
remains negative throughout the collapse indicating contraction of the collapsing 
configurations. 


5.2 Approach: II 

In this approach, we adopt an appropriate approximation procedure to examine 
the collapsing configurations with k /= 0. We set rs = 1 in equation Eq. (S3 
rewriting it as 

2(\/ (1 - k)£ + r])(y /(1 - k)£ + fj) - 3(y / (1 - k)£ + i)) 2 
-2 k(sj (1 - fc)£ + ii)i + 4£fc((l - k)i + \/(l- k)rj) 

-k(y/(l - k)£ + r,) 2 + 3fc 2 / 2 = 0. (51) 

Since a closed form solution of Eq. m is not available, we generate an approxi¬ 
mate solution by assuming 

V(t) = Pe nt , t(t) = eh(t), (52) 

with 0 < e << 1. By considering terms up to O(e), an approximate solution of 
the Eq. (1511) is obtained in the form 

e"V ( 1 ~k)(k + n 2 )P 


(53) 


L = y k{—4 + 5k) + 6 a/ (1 — k)kn — 5(—1 + fc)n 2 , 

and P, n, k, C, D are constants which should be fixed based on physical require¬ 
ments. 

The behaviour of the physical quantities of the collapsing configuration for this 
particular model has been examined for configurations with the initial mass and 
radius 5 Mq and 15 km, respectively. The evolution of the physical quantities for 
specific choices of parameters has been described graphically. Fig. [5] shows that 
the surface density p(rz, t ) starts from a finite value and increases as time evolves. 
The heat-flux q is zero initially and remains positive as time evolves (Fig. [ 6 ]). The 
total mass m(rs,t) and geometric radius ( r/A(r,t))s decrease from their initial 
values as time progresses as shown in Fig. [7] and [HJ respectively. 




2[fc 2 + n 2 + k(— 1 + n(y / (l — k) — n)\ 

(fc+3Ti\/(l-fc)-L)t —fe)+L)t 


+e 


Ce 


2 a /(1 -fc) 


+ De 


Vo-*) 


where 


6 Discussions 

In this work, we have described relativistic conformally flat solutions for spherically 
symmetric fluid configurations, undergoing dissipative collapse. It is not difficult 
to find a similarity between these solutions and the solutions obtained earlier by 
Herrera et al[ TJj which follow by adopting a different procedure. In the Herrera 
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et al\ T3] approach, conformal flatness is ensured by equating the Weyl tensor of 
the space-time to zero. The choice of the space-time in our approach is such that 
this condition is fulfilled on adopting a space-time with geometry conformal to 
a conformally flat space-time. The RW spacetime is conformally flat and further 
characterized by the feature that its spatial sections are homogeneous. Accordingly 
the space-time coordinates of the collapsing configuration are such that the spatial 
geometry is conformally homogeneous and continues to be so throughout collapse. 
This approach brings out the nature of impact of the spatial geometry determined 
by the curvature parameter k on the evolution of the collapse in the presence 
of dissipation. Further the approach also indicates the inter-dependence between 
the inhomogeneities in the matter density distribution and presence of dissipative 
processes during the evolution of the collapse. The mass of the particular collapsing 
configuration of case k = 0 of Approach I, tends to zero as collapse progresses 
indicating that all energy is radiated away as heat as in a case indicated earlier by 
Banerjee et aim ■ However, in the case discussed in Approach II, it is not so. 

To summarize, the class of conformally flat solutions reported in this article 
have been utilized to generate physically viable models of spherically symmetric 
gravitationally collapsing systems dissipating energy in the form of radial heat 
flux. As for the future investigations, since Maxwell’s equations are conformally 
invariant, it would be interesting to incorporate charge into our model and inves¬ 
tigate the impact of electromagnetic field on the overall evolution of the system. 
This issue will be taken up elsewhere. 
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Fig. 1 Evolution of surface density. (We have assumed k = 1, n = 10 and Q = 1.) 



Fig. 2 Evolution of heat-flux. (We have assumed k = 1, n = 10 and Q = 1.) 



Fig. 3 Evolution of mass. (We have assumed k = 1, n = 10 and Q = 1.) 
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Fig. 4 Rate of expansion. (We have assumed k = 1, n = 10 and Q = 1.) 



Fig. 5 Evolution of surface density. 
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Fig. 6 Evolution of heat-flux. 


m(r 2 ,t) 



Fig. 7 Evolution of mass. 
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Fig. 8 Evolution of proper radius. 






